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Manipulating the orbital state in a strongly correlated electron system is of fundamental and technological 
importance for exploring and developing novel electronic phases. Here, we report an unambiguous 
demonstration of orbital occupancy control between t 2g and e g multiplets in quasi-two-dimensional 
transition metal oxide superlattices (SLs) composed of a Mott insulator LaCo0 3 and a band insulator 
LaAl0 3 . As the LaCo0 3 sublayer thickness approaches its fundamental limit (i.e. one unit-cell-thick), the 
electronic state of the SLs changed from a Mott insulator, in which both f 2 ^and orbitals are partially filled, 
to a band insulator by completely filling (emptying) the t 2g (e g ) orbitals. We found the reduction of 
dimensionality has a profound effect on the electronic structure evolution, which is, whereas, insensitive to 
the epitaxial strain. The remarkable orbital controllability shown here offers a promising pathway for novel 
applications such as catalysis and photovoltaics, where the energy of d level is an essential parameter. 

Transition metal oxides (TMOs) offer a considerable number of emergent physical phenomena, such as 
superconductivity, magnetism, and (multi)ferroicity. These phenomena originate from the complex inter- 
play among charge, spin, orbital, and lattice degrees of freedom associated with the d states of transition 
metals 1 " 3 . Due to their inherently coupled nature, even a small change in the orbital degree of freedom can bring 
about drastic evolution of the underlying electronic structure 2 " 5 . Therefore, finding a tuning knob for controlling 
the relative occupancy of d orbitals has been a central issue to tailor the intriguing physical properties of TMOs. 
The outcome of such a study, thus, not only deepens our fundamental understandings but also allows us to 
envision novel TMO -based electronic devices 3 " 5 . 

Recent advances in atomic-scale synthesis of TMO heterostructures and artificial superlattices (SL) offer 
unprecedented opportunities for the control of orbital states that cannot be realized within bulk counterparts. 
Examples include (1) controlling the physical interaction and relevant orbital energy levels at the interface 
between a cuprate and a manganite, resulting in a systematic tuning of the orbital occupation 5,6 ; (2) inducing 
high-T c cuprate-like orbital states in a LaNi0 3 /LaAl0 3 perovskite heterostructure by controlling the electron 
occupation of particular d states with the e g symmetry (x 2 -y 2 and 3z 2 -r 2 ) 7 ~ 9 ; and (3) modifying the e g orbitals 
systematically to optimize the catalytic activities at the surface of TMOs 10 " 12 . The charge transfer needed for a 
catalytic reaction was found to be facilitated by the d-level near the chemical potential. 

Up to date, however, the orbital control has been rather limited to a subset of d electronic states, i.e. either 
within t 2g or e g states. Despite a huge number of reports on the control of the electron number within the 
respective t 2g or e g states (e.g. Refs. 7,9), the artificial control of relative electron populations across the two 
subsets has not been realized because the t 2g -e g separation, called lODq, is usually a few eV (in typical perovskite 
TMOs), much larger than the energy scale we can engineer e.g., with heterostructuring. In this Report, we 
demonstrate that a dramatic control over the occupation between the t 2g and e g orbital states can be achieved 
in quasi-2- dimensional (2D) artificial SLs of a TMO heterostructure. The SLs are composed of LaCo0 3 (LCO) as 
an active layer and LaAl0 3 (LAO) as a spacer layer. LCO is a unique perovskite which has a "twin" of quantum 
many-body states. In the case of bulk LCO, a low spin (LS) state (t 2g ; S = 0) is the ground state and a high spin 
(HS) state (t 2g 4 e g 2 ; S = 2) is the first excited state, with an energy difference less than 50 meV 13 " 15 (see 
Supplementary Information for details of the discussions on the intermediate spin state 16 ). At room temperature, 
bulk LCO has a mixed spin state of 40% HS + 60% LS due to thermal excitation 13 " 15 ; namely, ~1 electron (hole) 
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Figure 1 | Reducing dimensionality to control the c/-orbital population in 

(LaCo0 3 ) 1 j(LaA10 3 ) 11 . Two competing spin-orbital configurations, i.e. 
low-spin (S = 0; t 2g 6 ) and high-spin (S = 2; t 2g 4 e g 2 ), are shown for the Co 3+ 
ions. The two spin states in 3D bulk LaCo0 3 are nearly degenerate, 
resulting in a mixed-spin state at room temperature, whereas the low-spin 
state can be induced in quasi-2D superlattices by changing the energies of 
the 3z 2 -r 2 and xz/yz states. 

occupies the e g (t 2g ) orbital in bulk LCO. Therefore, the room tem- 
perature electronic state of bulk LCO has been considered as a Mott 
insulator 1417 . The energies of the twin states are in delicate balance 
and thus the electron populations of those states could be tuned with 
a minimal energy cost by applying strain or controlling dimension- 
ality through heterostructuring. 

Results 

Figure 1 shows a schematic diagram of our approach to control the 
inter-multiplet orbital occupation in LCO. At the ultrathin limit, one 
of the e g orbital levels, namely, 3z 2 -r 2 (z is the direction of the surface 
normal, i.e. along the crystallographic [001] direction) can be modu- 
lated to change the insulating nature of LCO. In particular, the 
quantum confinement effect can selectively increase the energy of 
orbitals concentrated along the z direction (3z 2 -r 2 and xz/yz). 
Moreover, in heterostructures with a spacer layer such as LAO, cer- 
tain chemical interactions at the interface between the LCO layer and 
the spacer layer can influence the strength of the Co 3d(3z 2 -r 2 )-0 2p 
hybridization 3 ' 818 . Such interactions can affect the total energy of the 
HS state while keeping the total energy of the LS state almost 
unchanged. Thus, heterostructuring as in a quasi-2D oxide system 
offers a means to modify the delicate population balance between two 
quantum many-body states (HS and LS). 

This scenario in ultrathin LCO is clearly supported by the results 
of the dynamical mean field theory (DMFT) calculations 19 . Figure 2 
comparatively shows the electronic structures of an ideal bulk (3- 
dimension; 3D) and a single layer (2D) LCO system (for details of the 
calculations, see Supplementary Information 16 ). The total/orbital- 
resolved density of states (DOS) for the 3D and 2D systems are 
presented in Figures 2a and 2b, respectively. The corresponding 
momentum-resolved spectral functions are shown in Figures 2c 
and 2d. The energy distribution of the orbital states shows broad 
Co e g (narrow Co t 2g ) bands located between +1 and +4 eV ( — 1 
and + 1 eV). O 2p bands are mostly located below —2 eV. In the case 
of 3D LCO, a portion of the e g orbital is occupied (the red-shaded 
area in the DOS near E = —0.8 eV), while a portion of the t 2g orbital 
is unoccupied (the blue-shaded area in DOS just above the Fermi 
level). This configuration manifests a mixed spin state in 3D LCO at 
room temperature 1415 . In contrast, in the 2D case, it is clearly shown 
that the e g bands are empty and the t 2g bands are almost fully occu- 
pied, suggesting an electron transfer from the e g orbital to the t 2g 
orbital. 

It should be noted that the total DOS near E = —0.5 eV is signifi- 
cantly enhanced as the dimensionality is reduced (3D — > 2D; from 
Figure 2a to Figure 2b). This is primarily due to the increase in the 
number of oxygen ions that participate in the O p - Co d hybridiza- 
tion in 2D LCO; in 3D LCO, the total DOS comprises 3{Op G ) + 6(0 
p n ) + 2(Co e g ) + 3(Co t 2g ), while in 2D LCO, 4(0 p a ) + S(0 p n ) + 




-3 -2 -1 0 1 2 3 4-4 -3 -2 -1 0 1 2 3 4 
Energy (eV) Energy (eV) 




Figure 2 | Dynamical mean field theory results. Total/orbital-resolved 
density of states (DOS) and spectral functions as a function of energy and 
momentum for (a, c) 3D and (b, d) 2D configurations, respectively. Each 
orbital-resolved DOS is normalized so that the total area becomes unity. A 
yellow arrow in (c) denotes the nearly dispersionless t 2g band. Possible 
optical transitions a, /?, and y are indicated by arrows. 

(x 2 -y 2 ) + (3z 2 -r 2 ) + (xy) + 2(yz/zx). Therefore, the charge transfer 
nature becomes much enhanced in 2D LCO even with a small O p n 
DOS at £ = -0.5 eV. 

The redistribution of the d-orbital occupation in dimension- con- 
trolled LCO yields a dramatic change in the nature of the electronic 
ground state from a Mott insulator to a band insulator. In the spectral 
function of the 3D LCO sample (Figure 2c), nearly non-dispersive 
and incoherent states (denoted by a yellow arrow) are observed just 
above the Fermi energy. These states have the Co t 2g orbital character, 
reflecting that the 3D LCO is a Mott insulator, as we denoted the 
optical transition across the Mott gap by the arrow a. In contrast, the 
2D LCO sample is found to be a band insulator, in which the gap at 
the Fermi level is defined by the energy separation between the fully 
occupied t 2g and empty e g bands (instead of t 2g ); the corresponding 
optical transition is represented by the arrow It is also shown that 
the lower e^band edge along the T-X direction shifts to higher energy 
in the 2D LCO case. The blueshift of the transition /? is presumably 
caused by the fact that the coordination of the Co 3+ ion with the 
upper and lower layers is absent, and the hopping in the out-of-plane 
direction is suppressed. The dimensional control can stimulate the 
electronic quantum phase transition as well as the sizable redistribu- 
tion of the orbital occupation in the 2D LCO sample. 

In order to experimentally validate the theoretical prediction on 
the control of orbital occupation, we designed TMO SLs composed of 
ultrathin LCO layers embedded within the large band-gap insulator 
LAO. Systematic dimensional control was conducted by reducing the 
periodicity of SLs below a few atomic unit-cells; namely, (LCO)„/ 
(LAO) n SLs (n = 2, 6, and 10). The LCO SLs were fabricated on 
NdGa0 3 (NGO) (110) substrates using pulsed laser epitaxy (see 
Supplementary Information 16 ). Figure 3 shows x-ray diffraction 
reciprocal space maps of the three samples (n = 2, 6, and 10) near 
the (103) reflection of the substrate. The indices on the satellite peaks 
reflect the periodicity of the SLs. For more detailed structural ana- 
lyses, see Supplementary Information 16 . 

It is clearly shown that all of the H values in the HiCL-coordinate of 
the SL peaks (red spots) were unity, indicating that the in-plane 
lattice constants were coherently maintained with respect to the 
substrate. This result suggests negligible influence of strain among 
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Figure 3 | X-ray diffraction reciprocal space maps of the three 
(LaCo0 3 ) n (LeLA\0 3 ) n superlattices with n = 2, 6, and 10, near the (103) 
reflection of the substrate. The indices on the satellite peaks reflect the 
superlattice periodicity. All of the H values in the ffKL-coordinate of the 
superlattice peaks (red spots) were unity, suggesting negligible influence of 
strains on the atomic structure of the LaCo0 3 superlattices. 

the LCO SLs. Note that the lattice mismatch of bulk LCO (pseudo- 
cubic) and LAO on a NGO substrate is 1.6% and 1.9%, respectively. 
However, it has also been shown that the electronic structure in LCO 
thin films is rather insensitive to the lattice strain, while the magnetic 
properties showed strong dependence 20 . The insensitiveness against 
the strain might be due to the peculiar electronic structure of LCO. As 
shown in Figure 2, there exists O 2p density of states near the Fermi 
level. According to our DMFT study on the strain dependence (not 
shown), the strain mostly influences the O 2p energies but not much 
the Co 3d energies. This suggests that valence, and consequently, 
ionic size of Co 3+ can hardly change under external stresses. 
Therefore, we can tell unambiguously that the evolution in the elec- 
tronic structure found in this work, is originated mainly from the 
dimensionality rather than the strain effects. 

The evolution in electronic structure predicted by the DMFT cal- 
culations is verified experimentally using optical spectroscopy. 
Figure 4 shows optical conductivity (0(00)) of the three SLs. For 
comparison, o(co) of a 13 nm-thick LCO film is attached. We 
observed strong absorption features at hco ~3 eV (thin gray lines) 
and 1.5 eV (thin yellow triangles), together with a weak absorption 
feature at hco ~ 0.5 eV (filled red triangles). Based on the results of 
the DMFT calculations (Figure 2), we can attribute these features to a 
O p-Co d charge transfer transition (y ~ 3 eV), a Co d(t 2g )-d(e g ) 
transition (/? ~ 1.5 eV), and a Co d(t 2g )-d(t 2 g) transition (a ~ 
0.5 eV), respectively. It is clearly observed that as n decreases, the 
weak absorption at to —0.5 eV (a) is systematically suppressed and 
disappears at (LCO) 2 /(LAO) 2 . Based on the DMFT results, the dis- 
appearance of the excitation a can be explained in terms of the cross- 
over of the spin-orbital ground state from a mixed HS + LS 
configuration to a LS configuration with decreasing the dimension, 
which is generically accompanied by a crossover from a Mott insu- 
lator to a band insulator; the optical transition between the Hubbard 
bands of Co t 2g states disappears with decreasing n, and the optical 
gap is defined by the transition from fully occupied Co t 2g bands to 
empty e g bands. Also, the blueshift of transition predicted by the 
DMFT calculation results, is indeed observed in the experimental 
0(00), as denoted by the yellow triangles in Figure 4. The excellent 
agreement between theoretical and experimental results confirms 
that our SLs indeed undergo the dimensional crossover of the elec- 
tronic structure, as illustrated in Figure 1. 

Discussion 

Polarization-dependent O iC-edge X-ray absorption spectroscopy 
(XAS) data further reveal the details of orbital occupation change 
in our (LCO) J (LAO) n SLs with variation in the dimensionality (see 
the Supplementary Information for details 16 ). Since the O iC-edge 
XAS probes the transition to the unoccupied O 2p-Co 3d hybridized 
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Figure 4 | Optical conductivity spectra of (LaCo0 3 ) 11 (LaA103) 11 
superlattices (n = 2, 6, and 10) and a 13 nm-thick LaCo0 3 film. With 
decreasing n y the optical transition within Co t 2g (i.e., kg—* h g ) was 
suppressed, opening a charge gap. Thin lines are fitting curves highlighting 
each peak. 



states, we can monitor the change in the 3d orbital occupancy of Co 3+ 
ions directly. Figure 5a shows the isotropic XAS spectra ([(£ // c) + 2 
X (E II ab)]!3) measured at room temperature. The spectrum from a 
13 nm-thick LCO film (gray dots) is included in the figure for com- 
parison. The main features near the photon energy hco = 529 eV and 
the shoulder near hco = 528 eV are attributed to the unoccupied Co 
e g and Co t 2g states, respectively 15 ' 21,22 . Features near hco = 532 eV 
reflect the Al 3sp states in the LAO substrate. In contrast to the case of 
thick LCO, the unoccupied Co t 2g state nearly disappears with 
decreasing n, while the main structure of the unoccupied Co e g state 
is clearly enhanced as shown by the arrow in Figure 5a. This spectral 
evolution provides concrete evidence of the electron transfer from 
the e g state to the t 2g state. 

Moreover, the polarization dependence in the XAS data reveals 
that the e g -to-t 2g electron transfer can be facilitated by an increase in 
the 3z 2 -r 2 orbital energy. Figure 5b highlights the anisotropy in the 
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Figure 5 | (a) O iC-edge XAS isotropic spectra of (LaCo0 3 ) n (LaA10 3 ) n 
superlattices (n = 2, 6, and 10) and the thick LaCo0 3 film taken in a 
fluorescence yield (FY) mode. The intensity of the t 2g (e g ) feature near 
528 eV (529 eV indicated by a vertical dotted line) decreases (increases) 
with decreasing n y suggesting a dramatic transition to a low-spin state. The 
measurement geometry is appended in the inset, (b) Polarization- 
dependent spectra and (c) the X-ray linear dichroism (XLD) spectra 
highlighting the orbital character. Features near 532 eV reflect the Al 3sp 
states in the LaA10 3 substrates. 
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electronic structure. The polarization- dependent XAS spectra for n 
= 2 and n = 10 SLs are shown in Figure 5b. The (E // c) data are 
shown as dots while the (E // ab) data are shown as solid lines. The e g 
peak in the (E II c) spectra are located at slightly higher energies 
compared with their respective peak in (E // ab) counterparts, sug- 
gesting that the state is at a higher level than the x 2 -y 2 state. 

The corresponding X-ray linear dichroism [XLD = (E // c)-(E II 
ab)] is displayed in Figure 5c with a dotted zero-line for guidance to 
eyes. The spectral difference is more clearly observed in the XLD 
spectra. The positive (negative) signs in the XLD indicate the dom- 
inance of the 3z 2 -r 2 (x 2 -y 2 ) character. The XLD signal is enhanced 
with decreasing n from 10 to 2. Such anisotropy in the e g orbitals 
manifests 'lifting' of the 3z 2 -r 2 level under low dimensionality 5,812,18 ' 23 . 
The energy difference between the 3z 2 -r 2 state and the x 2 -y 2 state is 
estimated to be ~1 eV, which is very close to the energy distance 
between the leading peak positions in the unoccupied ^-r 2 and x 2 -y 2 
DOS in the DMFT result (see Figure 2b). For more clear assignments, 
we also performed Co L 2 ,3-edge XAS and the corresponding XLD 
measurements. The analysis of the Co L 2y3 -edge data using cluster 
model calculations confirmed the XLD results of the changes in the e g 
orbital states (see Supplementary Information for more details 16 ). 

Our combined theoretical and experimental results clearly 
revealed how the d- orbital occupation with different symmetries 
(t 2g and e g ) can be manipulated by two-dimensional layering of 
atomically thin Mott insulator-band insulator SLs. By controlling 
the delicate competition between the crystal field splitting and elec- 
tron correlation effect, we found that a subtle change in the 3z 2 -r 2 
orbital could be induced as a consequence of a d-orbital reconstruc- 
tion from a t 2g ' 3 e g s state to a complete t 2g state. The lattice relaxation 
for all of the LCO SLs was not observed, implying that the evolution 
in the electronic structure is indeed driven by the reduced dimen- 
sionality rather than the lattice strain effect. The reduced dimension- 
ality altered the insulating nature of LCO with a spin state crossover 
from LS+HS to LS. Therefore, our demonstrated controllability of 
electron occupation between the t 2g and e g orbital states with atomic - 
scale SLs manifests that an accurate growth control of oxide hetero- 
structures can open up many opportunities for discovering new 
physical properties or functionalities such as catalysis and photovol- 
taics, where the energy of d level is an essential parameter. 

Methods 

Sample fabrication. (LCO)„/(LAO) M SLs (n = 2, 6, and 10) were grown on NdGa0 3 
(110) substrates by pulsed laser epitaxy at 700°C under an oxygen atmosphere 
(100 mTorr). A KrF excimer laser (/I = 248 nm) with a laser fluence of ~ 1 Jcm" 2 was 
used to ablate sintered LCO and single-crystalline LAO targets. The crystallinity of 
the fabricated SL samples were confirmed by lab x-ray diffraction (XRD) and hard 
XRD at the BL14B1 beamline of the Shanghai Synchrotron Radiation Facility in 
China. 

Optical spectroscopy. Optical conductivity spectra (a(co)) were obtained by 
reflectance measurements in the energy region of 0.1-1.0 eV and variable angle 
spectroscopic ellipsometry (JA Woollam V-VASE) in the energy region of 0.74- 
4.5 eV. The o~(co) of the LCO layers in the SLs were obtained by optical simulations 
implemented in the W- VASE software. The optical constants of bulk LAO were 
measured independently for the SL optical model calculations. 

X-ray absorption spectroscopy. Soft X-ray absorption spectroscopy (XAS) 
measurements were performed at the 2A beamline at the Pohang Light Source (PLS) 
in Korea, which is equipped with an elliptically polarized undulator (EPU). The 
energy resolution was approximately 0.1 eV. We measured the absorption 
coefficients in both the fluorescence yield (FY) and total electron yield (TEY) modes. 
Note that we have not observed any discernable difference between the TEY and FY 
data in the energy region of interest in the O K-edge XAS spectra. Here, we changed 
the polarization of light by tuning the phase and gap in the EPU while fixing the angle 
of incidence with respect to the surface normal to the SLs to 70°, rather than by 
rotating the samples. This process circumvents any issues associated with the finite 
probing depth of the soft X-rays. The spectra of E II a reflect only in-plane orbital 
hybridization, whereas the spectra of E II n reflect mostly (88%) perpendicular 
hybridization and partially (12%) in-plane orbital hybridization. The isotropic term 
([(£ llc) + 2X(E II ab)]/3) and the XLD ((£ // c) - (E II ab)) at the Co L 2;3 -edge were 
deduced from the raw data of (E 1 1 a) and (E I In) using the relationships of (E 1 1 a) = 
(E II ab) and (E II n) = sin 2 70° (E II c) + cos 2 70° (E II ab), respectively. The spectra 



were normalized maintaining the total spectral weight after subtracting the 
contributions of the NdGa0 3 substrates. 
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